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ABSTRACT 

Using Spitzer archival data from the SAGE (Surveying the Agents of a Galaxy's Evolution) pro- 
gram, we derive the Cepheid period-luminosity (P-L) relation at 3.6, 4.5, 5.8 and 8.0 microns for 
Large Magellanic Cloud (LMC) Cepheids. These P-L relations can be used, for example, in future 
extragalactic distance scale studies carried out with the James Webb Space Telescope. We also derive 
Cepheid period-color (P-C) relations in these bands and find that the slopes of the P-C relations are 
relatively flat. We test the nonlinearity of these P-L relations with the F statistical test, and find 
that the 3.6/xm, 4.5/xm and 5.8/im P-L relations are consistent with linearity. However the S.O/xm P-L 
relation presents possible but inconclusive evidence of nonlinearity. 
Subject headings: Cepheids — distance scale 



1. INTRODUCTION 

The Cepheid period-lmninosity (P-L) relation is 
an important component in extragalactic distance 
scale and cosmological studies. The most widely 
used P-L relations in the literature are obtained 
from Large Magellanic Cloud (LMC) Cephe ids in 



optical BVI (e.g., see 



Tanviij [199 9': 'Udalski et_J] 



Madore fc FreedmanI 
' 1999at ISandage et al.l 



1991 



2004 



Kanbur fc Ngeow 200 6) and ne ar infr a-red (NIR) Jgj^ 
(e.g., see [M adore fc FreedmanI [1 991: Gie ren et aLlll998t 
Groeriewegen. 2 000: ^N ikolaev et al.. .2004i : iPersson et ail 



2004 iNgeow et al.l 
( Bauer et al.l Il999f) 
"H 1120071 ^ 



20041 INgeow et al. 200^ bands. In addition, there are 
also some LMC P-L rela tions obt ained for MACHO 
VmachoRmacho (Nikolaev et all 
[20© and EROS VEmsMEm^ 
bands. Recently, INgeow fc Kanbud (|2007| ) applied a 
semiempirical approach to derive the LMC P-L rela- 
tion in Sloan ugriz bands: still within optical and/or 
NIR regimes. Therefore, Cepheid P-L relations are well- 
developed for wavelengths ranging from optical to NIR. 

In contrast, there are currently no P-L relations avail- 
able for wavelengths longer than the K band in the lit- 
erature. The main motivation for having a P-L relation 
at these longer wavelengths is in order to apply it in 
future extragalactic distance scale studies. The Near In- 
frared Camera (NIRCam) and the Mid- Infrared Instru- 
ment (MIRI), that are scheduled to be installed on the 
James Webb Space Telescope (JWST), will operate in 
the NIR and mid-infrared: a wavelength range of 0.6- 
5 microns and 5-27 microns, respectively^. It is possi- 
ble that extragalactic Cepheids will be discovered and/or 
(re-)observed (for those galaxies that were observed by 
the HST Ho Key Project) by the JWST. Hence, P-L re- 
lations at longer wavelengths are needed in order to use 
JWST to derive a Cepheid distance to these galaxies. In 
this Paper, we derive the LMC P-L relations at 3.6, 4.5, 
5.8 and 8.0 microns from the Spitzer archival data. As 

^ See the links given in |http : //www, stscl . edu/jMst/lnstruinents/| 



far as we are aware, this is the first time such a relation 
has been derived. 

In addition to the P-L relations, we can also derive 
period-color (P-C) relations and construct color-color 
plots and the color-magnitude diagrams (CMD) for LMC 
Cepheids in these Spitzer bands. In Section 2 we discuss 
our data selection. In Section 3 we present our analy- 
sis and results for the P-L relations. Section 4 we show 
the P-C relations, the CMD and the color-color plot for 
the Cepheids in our sample. Our conclusion is given in 
Section 5. Extinction is ignored in this Paper because it 
is expected to be negligible in the Spitzer 's IRA C bands 
(hereafter /i?^C band). 

2. DATA SELECTION 

The SAGE (Surveying th e Agents of a Galaxy's Evo- 
lution, [Meixner|eF^n[2QQl)^ is a program to survey the 
LMC using the Spitzer satellite. As a result it has de- 
tected about 4 million sources in the LMC with Spitzer 's 
IRAC instrument (with an angular resolution of ~ 2 
arcsecond). The data are publicly available via the 
IRSA's Gator Catalog Query^. To find LMC Cepheids 
in the SAGE database, we first obtained the right ascen- 
sion (RA) and dechnation (DEC) of the LMC Cepheids 
from the OGLE (Optical G ravitational Lensing Experi- 
ment, [UdnikLit!aiJ[i999B) database. We only obtained 
the Cepheids that were classified as fundamental mode 
Cepheids (labelled as FU) from the OGLE database, 
where the classification was mainly based on the Wi 
P-L relation and the Fourier decomposition technique 
(jUdalski et al.lll999b[ ). The initial list of 771 Cepheids 
was cro ss-correlated with a list of "good" Cepheids 
given in 'Kanb ur fc Ngeowl (|2006l ). The details regard- 
ing selection criteria used to remove "bad" Cepheids, 
including possible overt one Cepheids, can be found in 
iKanbur fc Ngeowl (|2006[ ) and will not be repeated here. 



See http : //sage . stscl . edu/ index .php' 
^ See http : //Irsa. Ipac . caltech. edu/ap pllcat Ions/Gator /"j 
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TABLE 1 

Search Results from the SAGE Database. 



TABLE 2 
P-L Relations in ffij4CBAND. 



SAGE Catalog 



SAGE Archive 







886 




912 


-^nomatch 




13 




4 


-'Vflnal 




724 




733 


<d> 




0.777 ±0.015 




0.773 ±0.014 




St. 


deviation= 0.393 


St. 


deviation= 0.374 


< ARA > 




0.575 ±0.020 




0.582 ±0.019 




St. 


deviation= 0.531 


St. 


deviation= 0.518 


< ADEC > 




-0.088 ± 0.014 




-0.092 ±0.013 




St. 


deviation= 0.372 


St. 


deviation= 0.349 



Note. — 
arcsecond. 



Units for < d >, < ARA > and < ADEC > arc in 



SAGE Catalog 



SAGE Archive 



distance or separation (arcsec) 

Fig. 1. — The distribution of the "distance", or the separation, 
between the Cepheid locations and the locations of the matched 
sources. 

This left 627 OGLE LMC Cepheids. To increase the 
sample size and to extend the period coverage to longer 
period (OGLE Cepheids truncated at P ~ 30 days due 
to CCD saturat i on) , w e added non-OGLE Cepheids from 
the iSebo et al. ' (2002) catalog. As in previous studies 
(for example, sec Udal ski et al.lfl999a: Kanbur fc Ngeow 
[20MlNgeow"eral.H2005[ ISandage et al.i2004l and the ref- 
erence therein), we apphed a period cut of P > 2.5 days 
to our sample to avoid inclusion of overtone Cepheids. 
The number of Cepheids left in our sample is 737. 

The IRSA's Gator Catalog Query provides two SAGE 
datasets: the SAGE Winter '07 IRAC Catalog (here- 
after SAGE Catalog) and the SAGE Winter '07 IRAC 
Archive (hereafter SAGE Archive). The difference be- 
tween the SAGE Catalog and the SAGE Archive can be 
found in Meixner et al. (2006) and in the SAGE docu- 



ment , and will not be repeated here. The search on 
the SAGE Catalog and Archive were performed with a 
search radius of 3.5 arcseconds. The results are sum- 
marized in Table [TJ In this table, iVtotai is the total 
number of matched sources. However there are several 
Cepheids which do not have any matched sources within 
the search radius: the number of these Cepheids is given 

* See http : //Irsa. ipac . caltech. edu/ application s/Gator/] 
GatorAid/SAGE/SAGE_SSCdatadocument_delivered.pdf 



Band 


Slope 


Zero-Point 




N 




SAGE Catalog 


3.6/im 
4.5/im 
5.8/im 
8.0/im 


-3.265 ±0.017 
-3.211 ±0.017 
-3.158 ±0.027 
-3.031 ±0.048 


15.947 ±0.012 
15.922 ±0.012 
15.840 ±0.022 
15.609 ± 0.049 


0.104 
0.104 
0.169 
0.183 


613 
618 
534 
215 


SAGE Archive 


3.6/im 
4.5/im 
5.8/im 
S.Ofim 


-3.263 ±0.016 
-3.221 ±0.017 
-3.173 ±0.028 
-3.091 ±0.039 


15.945 ±0.012 
15.927 ±0.012 
15.850 ±0.022 
15.684 ± 0.036 


0.104 
0.103 
0.175 
0.193 


628 
635 
561 
319 



Note. — rr is the dispersion of the P-L relation. 

as A^nomatch- For thosc Cepheids with multiple matched 
sources, the sources with a minimum "distance" or sepa- 
ration from the input RA and DEC of the LMC Cepheids 
were selected to remove the duplicated SAGE sources. 
The number of matched sources left in the sample is 
summarized as Affinal in Table [1] The mean "distance" , 
< d >, between the input Cepheid locations and the 
locations of the matched sources, and the mean differ- 
ence in RA and DEC, ARA and ADEC respectively, 
for the remaining Cepheids are presented in Table [T] as 
well. Figure [1] shows the distribution of the "distance" or 
the separation between the remaining Cepheids and the 
matched sources from both of the SAGE Catalog and the 
SAGE Archive. 

Running the same query on SAGE Winter '07 MIPS 
24/im Catalog returned 12 and 25 matched sources with 5 
and 8 arcsecond search radius, respectively We therefore 
do not consider the detections from the MIPS catalog in 
this Paper. 

3. THE PERIOD-LUMINOSITY RELATION 

Initial plots of the P-L relations for all matched sources 
display a tight P-L relation with some obvious outliers. 
These outliers are probably due to the mis-match of the 
Cepheids with the input SAGE database, blending of 
other sources along the line-of-sight, or other physical 
reasons. These reasons are difficult to track down be- 
cause we only have the data from the publicly available 
database and not the source images. Therefore we apply 
an iterative outl ier removal algorithm (the sigma-clipping 
algorithm, as in lUdalski et al.llT999af ) to remove the out- 
liers. For each iteration, P-L relations are fitted to the 
data and outliers located more than 2.5cr away from the 
fitted regression lines are removed, where a is the dis- 
persion of the regression lines. This process is repeated 
several times until the solutions from the regression be- 
come stable. The rejected outliers are represented as 
open squares in Figure [2l 

Figure [2] displays the P-L relations in the ZR^C bands 
and Table [2] presents the results from fitted regression 
lines after outliers have been removed. When fitting 
the P-L relations, we do not stipulate constraints that 
the number of Cepheids should be the same in all four 
bands and/or a given Cepheid is detected in all four 
bands because the number of Cepheids in the 8.0/im 
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Fig. 2. — The P-L relations from the matched sources in the SAGE Catalog (left panel) and the SAGE Archive (right panel), after an 
iterative process to remove the outliers. The open squares and solid circles are for the rejected outliers and the remaining data points, 
respectively. The lines are the fitted P-L relations to the data as given in Table [2] 



band is much smaller than in other bands. Table [5] finds 
that the /i?^C band P-L relations obtained using data 
from the SAGE Catalog and the SAGE Archive are in 
very good agreement. A small discrepancy is seen for 
the 8.0/im P-L relation, perhaps due both to the differ- 
ent numbers of Cepheids in the two catalogs (there are 
^ 48% more Cepheids from the SAGE Archive than in 
the SAGE Catalog) , and also the lack of Cepheids fainter 
than m = 14mag. in the SAGE Catalog (see lower-left 
panel of Figure [2]). In Figure [31 we give the reported 
Icr errors on the magnitudes as a function of magnitude 
from the SAGE database. This figure clearly shows the 
truncation of m = 14mag. at 8.0/im for our Cepheid 
sample in the SAGE Catalog. Furthermore, there is a 
clear cut-off of the Icr errors for the 8.0/im band data^ 
from the SAGE Catalog, which is absent in other panels 
in Figure [31 The lack of ms.o > 14mag. and the cut-off 
of the Icr errors at 8.0/im band could be due to the more 
stringent crit eria for the SAGE C atalog than the SAGE 
Archive (see iMeixner et all [20061 and the SAGE docu- 
ment given in footnote [4]) . For these reasons, we only 
consider the data from SAGE Archive in the rest of this 
Paper. 

Figure [H presents the residual plots from the fitted 
P-L relation in each bands. The residuals for 3.6/im, 
4.5/ini and 5.8/im band are more or less evenly dis- 
tributed around the regression lines. However the resid- 
uals for 8.0/im band exhibit a deficit on one side of the 
regression at log(P) ^ 1.0. This could suggest that the 
8.0/im band P-L relation may not be linear, and this will 
be investigated further in Section [3.41 Interestingly, the 

^ This cut-off is still visible if we include a large number of non- 
Cepheid data from the SAGE Catalog. 



smallest dispersion of the P-L relations seems to occur 
at log(P) 1.0 in all four bands. 

3.1. Sensitivity of the P-L Relations with Various 
Period Cuts 

Even thought period cut of log(P) ~ 0.4 has been ap- 
plied to our sample to avoid the contamination of over- 
tone Cepheids at the short period end (see Section 2), the 
short period Cepheids in our sample may still be influ- 
enced by the overtone Cepheids, the increasing sensitiv- 
ity of blending and the increasing of measurement errors 
in the IRA C band (as shown in Figure [3]) for the faint 
(hence short period) Cepheids. These may bias the P-L 
relations when comparing to those presented in Table [2l 
and a cut at a longer period may be required to avoid 
the bias due to these effects. In Figure [51 we present the 
slopes and zero-points of the fitted P-L relations with 
various period cuts from log(P) = 0.4 to log(P) — 0.9 
to our sample. This figure suggests that the slopes are 
consistent (within the Icr error bars) at different period 
cuts up to log(P) ~ 0.65. A similar result is also found 
for the zero-points. This suggests that increasing the 
period cuts up to log(P) ~ 0.65 will not greatly affect 
the results presented in Table [2l and the various sources 
that could bias the P-L relations may not be important. 
However at log(P) ~ 0.65, the number of Cepheids in the 
sample decreases to ~ 50% in 3.6/im, 4.5/im and 5.8/tm 
bands, and to ~ 23% in 8.0/im band. Hence we retain 
our results given in Tabled in this Paper. 

3.2. Random Phase Correction and the P-L Relations 
in the IRAC Band 

In optical bands accurate mean magnitudes from 
the light curves are used to fit for the P-L relations. 
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Fig. 3. — Plots of the la errors on the magnitudes as reported in the SAGE database versus the magnitudes. The open squares and sohd 
circles are for the rejected outliers and the remaining data points, respectively. 
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Fig. 4. — The residuals from the P-L relations given in Table [2] 
(SAGE Archive) as a function of period. 

In contrast, the magnitudes obtained from the SAGE 
database are flux av eraged over several observations (see 
iMeixner et al.ll2006L for the details of the SAGE observ- 
ing strategy). Since information regarding the time of 
observation of each data point is not available to us, 
these averaged magnitudes may not necessarily corre- 
spond to the mean magnitudes of the Cepheids. Nev- 
ertheless, it is well-known that the amplitude of Cepheid 
light curves decreases as wavelength increases. From Fig- 
ure [6l it can be seen that the amplitude decreases from 
B to J band, and flattens out at H and K band. Ex- 
trapolating this trend to longer wavelengths may imply 
that the amplitudes in the ZR^C bands could either be 
the same as in HK band, or become smaller or larger as 
wavelength increases. Therefore, we assume that the am- 
plitudes will be small and these "random-phased" mag- 
nitudes arc close to the mean magnitudes in those bands. 
In JHK bands, there are well-developed methods to 
estimate the mean magnitudes from single epoch (or a 



few epochs) random phas e observations, a s prese nted in 
iNikolaev et "all (|2004D and lSoszvhski et all (|2005[ ). How- 
ever, both methods either require the epoch of obser- 
vation to be known or the existence of a template con- 
structed from well-observed light curves. Clearly, both 
requirements are not present for the data we studied in 
this Paper. Nevertheless, the use of random-phase JHK 
magnitudes from single epoch observations to derive the 
P-L relat ion has a precedent in the literature (see, for 
example, iGroenew egen 2000j . 

The influence of (single epoch) random phase magni- 
tudes is expected to be minimal in the K band. In order 
to test the validity of this assumption we use 2MASS K 
band data that are available from the SAGE database. 
We repeat the fitting of the K band P-L relation in the 
same way as for the IRAC bands. The results are pre- 
sented in Table [3l We compare our K band P-L rela- 
ti ons with the K band LMC P-L relation that is given 
in iFouque et al.l ()2007l ). There are two main reason s for 
selecting the P-L relation from iFouaue et all (120071) for 
our co mparison. Firstly, Cepheids used in iFouque et al.l 
()2007l ) are entirely based on the OGLE sample. This 
is similar to our sample that consists mostly of OGLE 
Cepheids. Secondly, random phase c orrection with the 
method described in iSoszvhski et aP (2005) and extinc- 
tion correcti on have been applied to the 2MASS K band 
data used in lFouque et al.l (|2007f ). This is in contrast to 
the data used in this Paper. 

Table [3] reveals that without applying any random 
phase correction, the slope of the K band P-L relation 
from the SAGE data base is identi c al to the K band 
P-L relation given in IFouque et al.l (|2007| ). This sug- 
gests that the random phase correction may not have 
a large influence on the K band slope, probably due 
to the small am plitude in the K b and. The differ- 
ence between the IFouque et al.l (|2007i) P-L relation zero 
points and the results found with the SAGE database is 
~ — 0.06mag. Assuming the mean extin ction toward the 
LMC is ~ O.lOmag. (for example, see iFreedman et "aD 
l200lh and using the t otal-to-selectiye ext inction (i?) of 
0.38 in the K band (jFouque et al.ll2007D . then about 
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Fig. 5. — The slopes (left panel) and the zero-points (bottom panel) of the fitted P-L relations as a function of adopted period cuts 
(log[P,„t]). 



TABLE 3 

Comparison of the K Band P-L Relation. 





Source 


Slope 


Zero-Point cr 


N 




Catalog 


-3.231 ±0.021 


16.050 ±0.016 0.140 


634 




Archive 


-3.229 ±0.021 


16.048 ±0.016 0.141 


642 


I 


F07 


-3.228 ±0.028 


15.989 ±0.006 0.136 


529 


Note. 


— (T is the 


dispersion of the P-L 


rela- 



H K 



X(//m) 

Fig. 6. — The typical ratio of the light curve amplitude (A) as a 
function of wavelength, using the V band amplitude as a reference 
(hence the ratio for V band is 1). The amplitude ratios a re from 
IFreedmanI ildM) for the B band. ITan vir (1997) & Ngeo w et aT] 
II2003I) for the / band, and ISoszvhski et al.i (|2005) for th e JHK 
band. There are two points for JHK band because Soszvfis ki et al.l 
II2005I) separated out the amplitude ratios at two period ranges. 

— 0.04mag. of the — O.OGmag. difference can be explained 
as due to the extinction. Therefore, we beheve that 
the IRAC band P-L relations presented here will not 
be strongly influenced by the lack of random phase or 
extinction corrections. 

3.3. Comparison of P-L Relations at Different Bands 

It has been well documented in the literature that 
the slopes of the P-L relation become progressively 
steeper from B to K band, while at the same time 
the dispersion of the P-L relation decreases (see, f or ex- 
ample, Madorc & Frccdnian 1991; Bcrdnikov ci^a]| |1996l : 
iCaputo et al.l i2000: Fi orentino et'al.ll2002t [2007f) . This 
is expected in part due to the black-body curves with 
Cepheid-like temperatures. From L oc R^T'^, the tem- 



tion. Catalog— SAGE Catalog 
FOT MFougue et al.l I I2007I) . 



Arehivo=SAGE Arehive; 



perature variation will dominate the luminosity varia- 
tion in the optical bands (see, for example. lCox|[T980f) 
and extends to JH band or even the K band. In con- 
trast, at K band and/or the wavelengths longer than the 
K band, the radius variation will dominate the tempera- 
ture variation. Since the period-radius (P-R) relation is 
independent of wavelength, the slope of the P-L relation 
is expected to reach a maximum value at some character- 
istic wavelength and remain constant as the wavelength 
increases. Similarly, the dispersion is expected to reach 
a minimum at the same characteristic wavelength and 
remain steady at longer wavelengths. In Figure [7l we 
compare the slope, the zero-point and the dispersion of 
the P-L relation at various bands. For illustration pur- 
poses, we adop t the empirical BVI a nd J HK band P-L 
relatio ns from ISandage et al.l (|2004l ) and iPersson et al.l 
(|2004l ). respectively. To extend to longer wavelengths, 
we also add P-L relations in the /i?^C bands from Table 
[2] (with results from SAGE Archive) in this figure. 

Figure [7| attests to the fact that the slope of the 
P-L relation is the steepest around the K to 4.5/xm 
band region, and becomes shallower at longer wave- 
lengths. In contrast, the zero-point of the P-L rela- 
tion is a monotonic function of wavelength. The dis- 
persion of the P-L relation displays a similar trend to 
the slope: the dispersion reaches a minimum around 
the K and 4.5/zm bands, and subsequently increases for 
longer wavelengths. A polynomial function in the form 
of Y{X) — ao + fti X [logjQ(A)]* was used to fit the 

data points in Figure [T] Here Y represents either the P- 
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Fig. 7. — Comparison of the slopes (top panel), the zero-points 
(ZP, middle panel) and the dispersions (bottom panel) for the em- 
pirical LMC P-L at different bands. The IRAC band data points 
are from the SAGE Archive results as given in Table[2] The dashed 
curves are the polynomial fits to the data points. 

L slopes, zero-points or dispersion. From the polynomial 
fits, the steepest P-L slope of —3.266 seems to be located 
at a wavelength of ^ 2.7 microns. The minimum disper- 
sion of ^ 0.109 occurs at the same wavelength. This 
suggests that the P-L relations from K band to 3.6^m 
band will provide a more accurate distance scale mea- 
surement than the optical bands. However, observations 
from optical bands are still required to detect Cepheid 
variables. Furthermore, the bottom panel of Figure [7] 
implies that the dispersions at 3.6/im and 4.5/im bands 
are similar to the JHK bands, suggesting also that omis- 
sion of random phase corrections may not be too serious 
an oversight. 

Figure [7] provides evidence that the slopes and disper- 
sions for the 3.6/im and 4.5/im band P-L relations are 
consistent both with each other and with the K band 
P-L relation. However, the slopes and dispersions for 
the 5.8/im and 8.0//m band P-L relations are shallower 
and larger, respectively than the theoretical expectation 
outlined previously. There are two possible causes for 
this. The first possibility is that measurement errors be- 
come larger toward the faint end of the 5.8/im and 8.0/im 
band P-L relations, as suggested in Figure [S] The second 
possibility is that the number of Cepheids in these two 
bands is less than those in the 3.6/im and 4.5/im bands. 
This is especially true for the 8.G/im band. Figure [8] 
presents the period distribution for our Cepheid sam- 
ple: fewer Cepheids are indeed detected for log(P) ^ 0.8 
in the 8.0/im band. This is also reflected in Table [5] 
where the number of Cepheids in the 8.0/im band (from 
the SAGE Archive), is about 43% to 50% less than 
the nu mber of Cepheids in other bands. iMeixner et al.l 
(|2006f ) reported the limiting magnitude in each band as 




Fig. 8. — The period distribution for the SAGE Archive samples 
that used to derive the P-L relation given in Table [2] 

'7^3X*4.5, 5.8, 8.0 = {18.348, 17.474, 15.164, 14.226}mag. 
This may result in fewer detections around m ~ 14mag. 
in the 8.0/im band. In the case of the 5.8/im band, Figure 
[8]also suggests that, for log(P) < 0.5, fewer Cepheids are 
detected when compared to the 3.6/im and 4.5/im bands. 

3.4. Nonlinear ity of the P-L Relations 

A number of recent studies have strongly sug- 
gested the LMC P-L relation in optical bands is 
nonlinear, in the sense that the relation can be bro- 
ken int o two P-L relation s separ at ed at /around 10 
days_ (iTammann fc Reindll l2002t iKanbur fc Ngeowl 
I2q04t iKa nbur fc Nge owl 12006 : 
iNgeow fc Kanbml l2006al lb 



12004'; 'Sa ndage et al 
Ngeow et al.| |2005 
Kanbur et al. 2007b: Ngeo w et all 



Ngeow et al.. (,2005.) and .Ngeow et al 
the JH band LMC P-L relations 
the K band LMC P-L relation 



2008 



In NIR, 
2008f l found that 
are nonlinear but 
is marginally linear. 
INgeow fc Kanbid (l2006br i outline a black-body argu- 
ment for their result that the P-L relation could be 
linear in K band but not in the optical and/or JH 
band. One possibility is that the temperature variation 
in Cepheid atmospheres, modulated at certain phases, 
periods and metallicities by the stellar photosphere- 
hydrogen ionization front i n teraction (IKanbur et al.l 
[200l IKanbur fc Ngeowl [20061 : IKanbur et al.ll2007al ). is 
responsible for the observed nonlinear P-L relation. 
Since the temperature variation for a black-body with 
Cepheid-like temperatures is minimal or even negligible 
at longer wavelengths, it is expected that the P-L 
relation becomes linear for wavelengths longer than the 
K band. Therefore, another main motivation for this 
Paper is to study the linearity /nonlinearity of the LMC 
P-L relations in the /i?j4Cband. 

To test the nonlinearity of the P-L r elations, we apply 
the F - test a s in our previous stu dies (jKanbur fc Ngeowl 
[200l [2006t INgeow et al.l [2005h . The detailed de- 
scri ption and forma lis m of the f-test c a n be found 
in [Wsberl (ll980D . IKanbur fc Ngeowl (|2004[ ) and 
INgeow et al.l (|2005l ). and will not be repeated here. Sim- 
ply speaking, in our f-test, the null hypothesis is that 
the data can be fitted with a single regression line, and 
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TABLE 4 

F-Test Results of the P-L Relations in IRAC Band. 



P < 10 days P > 10 days 

Band Slopes Zero-Pointy as Ns Slope/, Zero-Pointi, Nj^ F p{F) 

3.6/im -3.309 ±0.030 15.971 ± 0.019 0.101 550 -3.287 ± 0.079 15.986 ± 0.099 0.126 78 1.78 0.170 

4.5/im -3.255 ±0.029 15.946 ± 0.019 0.100 558 -3.331 ± 0.079 16.076 ± 0.098 0.122 77 2.81 0.061 

5.8/im -3.121 ±0.058 15.817 ± 0.039 0.176 471 -3.212 ± 0.093 15.893 ±0.118 0.175 90 0.55 0.579 

S.O^tm -2.858 ±0.094 15.517 ± 0.069 0.191 232 -3.364 ±0.114 16.018 ±0.143 0.186 87 5.90 0.003 



Note. — cr is the dispersion of the P-L relation. 

the alternate hypothesis is that two regression hnes sep- 
arated at 10 days are needed to fit the data. In our test, 
we set p{F), the probabihty of the observed F value un- 
der the null hypothesis, to be 0.05 (equivalently at the 
95% confidence level). This corresponds to ~ 3 for 
our data. Hence, F > 3 indicates that the null hyopthe- 
sis can be rejected at the 95% confidence level or more 
for the P-L relation under scrutiny. 

In Table [H we present the results from the F-test for 
the P-L Relations in ZR^Cband. The F-test results in- 
dicate that the P-L relation is linear in the 3.6/im, 4.5/im 
and 5.8/im band but not in the 8.0/im band. The linear- 
ity of the P-L relation in 3.6/im, 4.5/zm and 5.8/im band 
is expected from the black-body argument as outlined 
previously. However it is important to point out that the 
linearity /nonlinearity of the P-L relation in these bands 
does not necessary imply the P-L relation in the optical 
bands will be linear/nonlinear. The apparently nonlin- 
ear S.O/um P-L relation is puzzling. Removing the longest 
period Cepheid in the 8.0/im band sample still leaves a 
nonlinear result with a F value of 5.04. The relatively 
small number of Cepheids in 8.0/im band (see Figure [8]) 
may cause the apparent nonlinear result though the F 
test is sensitive to this. The similar trends of the la 
measurement error plot, as shown in Figure [3l and the 
similar dispersions for the P-L relations in the 5.8/im and 
8.0/im band, suggest the lack of short period Cepheids 
in 8.0/im band may be the reason for the nonlinear re- 
sult, because the F-test result finds that the 5.8/im P-L 
relation is linear. It is still inconclusive if the 8.0/im P- 
L relation is truely nonlinear or not and more data are 
needed in the future work to solve this problem. 

4. THE PERIOD-COLOR RELATION, THE COLOR-COLOR 
PLOT AND THE CMD 

In addition to the P-L relations derived in previous 
sections, the SAGE Archive also permits the derivation 
of P-C relations, color-color plots and the CMD. The 
resulting P-C relations from the data are presented in 
Figure [9] and Table [5l while the color-color plots and 
CMD are presented in Figures [10] and [11] respectively. 

The IRAC band P-C relations are found to be rela- 
tively flat as compared to the P-C relations in the op- 
tical band, with the mean color being close to zero (es- 
pecially for the TO4.5 — ms.g P-C relation). This is not a 
surprise given that the slopes and the zero-points of the 
P-L relations in these bands are similar (see Table [2] and 
Figure [7|) . In fact black-body curves with Cepheid- like 
temperatures predict the P-C relation should vanish at 
these bands. Table [S] also finds that some P-C relations 
are identical to each other, including the ma.g — ms.o & 
nT-s.s — rns.o pair, and the ma.g — TO4.5, 777,3. g — ™5.8 & 
1^5.8 ~ TTig.o P-C relations. Furthermore, all P-C rela- 
tions presented in Figure [9] display a very tight sequence 



TABLE 5 
P-C Relations in IRAC Band. 



Color 


Slope 


Zero-Point 


a 


N 


"13.6 


- m4.5 


-0.044 ±0.010 


0.017 ±0.007 


0.061 


595 


"13.6 


- Tns.S 


-0.048 ± 0.024 


0.050 ±0.019 


0.137 


504 


ms.e 


- rns.o 


-0.108 ±0.038 


0.190 ±0.034 


0.171 


287 


m4.5 


- ms.s 


-0.008 ±0.025 


0.039 ± 0.019 


0.140 


512 


m.4.5 


- ms.o 


-0.057 ±0.038 


0.169 ± 0.034 


0.170 


291 




~ ms.o 


-0.104 ±0.038 


0.186 ± 0.035 


0.181 


298 



Note. — cr is the dispersion of the P-C relation. 

for Cepheids with log(P) ^ 1.0. This is also seen in the 
CMD. In Table [5] we present the P-C relations separated 
at 10 days. The flatness of the P-C slope, the mean color 
of zero and the small dispersion of the P-C relation are 
clearly evident from this table for Cepheids with period 
longer than 10 days. 

At the short period end of the P-C relations, the dis- 
persion of the P-C relations gets broader (except for the 
fT^a.e ~ "^4.5 P-C relation) as period decreases. This fea- 
ture is also seen from Figure [TT] toward the faint end of 
the CMD. Probably this could be due to the relatively 
large measurement errors at the faint (or short period) 
end, because the measurement errors in color can reach 
up to ~ 0.3mag. as suggested from Figure |31 In addi- 
tion. Figures [9l and [TT] implv a lack of detections for the 
matched sources near log(P) ~ 0.5, especially for the 
P-C relations that include the 8.0/im band. This could 
cause the short period P-C relations to deviate from flat- 
ness and mean zero color as given in Table [5] Because 
of these reasons, we did not test the nonlinearity of the 
P-C relations with the F-test. 

Nevertheless, the CMD presented in Figure [TT] finds 
that the instability strip for LMC Cepheids is well de- 
fined in the /i?^C band, especially with the 7773.5 ~ "74.5 
color. The tightness of the 7773. g — TO4.5 color is also re- 
flected in Figure [TU] where the spread out of 7775.8 ~ ttis.o 
color is mainly from the short period Cepheids. The well 
occupied regions of the Cepheids in Figure [TO] and the 
well defined CMD suggest that Figure [TO] and [TT] can be 
used to identify Cepheids in future studies. 

5. CONCLUSION 

In this Paper, we derive P-L relations for LMC 
Cepheids in IRAC 3.6, 4.5, 5.8 and 8.0 microns bands. 
These P-L relations can be potentially applied to future 
extragalactic distance scale studies with, for example, the 
JWST. The data are taken from the Spitzer 's archival 
database from the SAGE program. After properly re- 
moving the outliers, the fitted P-L relations are presented 
in Table[2] We have tested the P-L relations with various 
period cuts and found that our results are insensitive to 
period cuts up to log Pent ^ 0.65. We also argue that 



8 



Ngeow & Kanbur 



I I I 



0.5 fr '■ .... 



T'^^^l^r^!^''^'^^^'-' ' 



I I I 



log(P) 



log(P) 



Fig. 9. — The P-C relations from the matched sources in the SAGE Archive after the removal of the outliers. 



TABLE 6 

The P-C Relations in IRAC Band Separated at 10 days. 







P < 10 days 








P > 10 days 






Band 


Slopes 


Zero-Point5 




Ns 


Slopci, 


Zero- Point i 




Nl 


"1.3.6 


- m.4.5 


-0.017 ±0.019 


0.001 ±0.012 


0.061 


521 


-0.031 ± 0.039 


-0.007 ±0.048 


0.059 


74 




- ms.s 


-0.117 ±0.050 


0.093 ±0.034 


0.146 


429 


-0.012 ± 0.037 


0.018 ± 0.046 


0.058 


75 


ms.e 


- ms.o 


-0.370 ± 0.097 


0.376 ± 0.072 


0.189 


213 


0.057 ±0.051 


0.001 ± 0.064 


0.079 


74 


m4.5 


- ms.S 


-0.118 ±0.051 


0.107 ±0.034 


0.149 


438 


0.027 ± 0.035 


0.015 ± 0.043 


0.053 


74 


"14.5 


- "18.0 


-0.335 ± 0.094 


0.364 ± 0.070 


0.185 


217 


0.096 ± 0.060 


-0.000 ± 0.075 


0.092 


74 


"16.8 


-ms.o 


-0.377 ±0.106 


0.381 ± 0.078 


0.204 


215 


0.075 ± 0.047 


-0.024 ± 0.059 


0.076 


83 



Note. 



- a is the dispersion of the P-C relation. 



Fig. 10. — The color-color plot from the 261 matched sources in 
the SAGE Archive, after the removal of the outliers, with detec- 
tions in all four bands. 

the random phase corrections may not be important for 
IRAC band P-L relations. When comparing P-L rela- 
tions from B to 8.0/im bands, the slope of the P-L rela- 
tion appears to be the steepest around K band to 3.6/im 
band, while the dispersion of the P-L relation reaches a 
minimum between those bands. The shallower slopes and 



larger P-L dispersions in the 5.8/um and 8.0/xm band are 
in contrast to the theoretical expectation. This could be 
due to the smaller number of Cepheids and larger mea- 
surement errors toward the faint end in these two bands. 
We also test the nonlinearity of the P-L relations in the 
IRAC band using the F statistical test. As expected, 
the F-test results show that the P-L relations are lin- 
ear in 3.6, 4.5 and 5.8 microns band, but the 8.0/um P-L 
relation is found to be nonlinear. However, the nature 
of the nonlinear 8.0/im P-L relation is still inconclusive. 
For the P-C relations, it was found that the slopes of the 
P-C relation are relatively flat in the /i?^C bands. Fi- 
nally, the LMC Cepheids show a well-defined instability 
strip in the CMD and clustered in a small region in the 
color-color plot. This can potentially be used to iden- 
tify Cepheids observed in the ZR^C bands. Even though 
there may be some associated problems for the 8.0/Ltm 
P-L relation, the 3.6/tim, 4.5/im and perhaps the 5.8/xm 
P-L relations can still be used in future distance scale 
studies. 
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